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Abstract 

Some of the progress in determining the phase boundaries of the 
nuclear phase diagram, the location of the critical point of the nuclear 
fragmentation phase transition, and the values of the critical exponents 
of this transition is reviewed. The recently postulated connection of 
the size of the largest fragments as a function of their rank with Zipf’s 
law known from linguistics is discussed, and an exact formula is de¬ 
rived. We show that Zipf’s law is not strictly observed in nuclear 
fragmentation, but that one can derive useful approximations. 

One of the premier goals of nuclear physics during the last two decades 
has been the determination of the phase diagram of nuclear matter and finite 
nuclei. From basic many-body physics, from symmetry arguments, from 
general considerations of the nucleon-nucleon force and of QCD, as well as 
from the analysis of experimental data, we have a fairly well established idea 
what the gross features of this phase diagram are. There are two distinct 
phase transitions, one between hadrons and quark/gluons, and the other 
between a nuclear “liquid” and a vapor of hadrons, see fig. ^ 

In order to learn about the parameters of the nuclear equation of state, 
we have to overcome several severe challenges. First and foremost, we cannot 
prepare our system at a given value of the state variables and keep it there 
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for an extended time interval. The only way in which we can excite our 
nuclear system is by colliding finite nuclei with each other or with small 
probes such as individual protons or pions. 
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Figure 1: Left: General features of the nuclear matter phase diagram. Right: 
circles: experimental data for the energy spectra of fragments produced in 300 
GeV proton-induced fragmentation events of Xe histograms: calculations with 
different values of the freeze-out density. The three cases on the right correspond 
to the areas of the phase diagram indicated on the left. 


A lot of effort has been devoted to constructing ways in which to provide 
reliable and necessarily indirect measurements of the state variables. A pri¬ 
ori, it is an open question whether thermodynamic equilibrium is achieved, 
and if equilibrium is achieved, at what time during the course of the reaction 
this happens first, and at what time different degrees of freedom freeze out. 
Further, it is unclear on which path the excited nuclear system migrates 
through the above phase diagram. 

However, the last few years of investigation have given rise to the hope 
that we can use nuclear reactions at certain energies, impact parameters, 
and target/projectile combinations to enable the excited nuclear system to 
pass sufficiently close to the critical point and stay there long enough for 
the final fragmentation pattern to contain evidence on the location of the 
critical point in the p/T diagram and perhaps the universality class of the 
phase transition. 

But how can a system with numbers of constituents of only several hun¬ 
dred at most can enable us to make reliable statements about a phase tran¬ 
sition that is defined strictly only in the infinite size limit? We would like 








to argue that nuclear physics has a privileged position and can provide ex¬ 
perimental and theoretical answers to this question that are relevant for 
many-body science in general. 

In order to probe the fragmentation phase transition one can either use 
symmetric heavy ion collisions at beam energies per nucleon around the 
Fermi energy or by using protons or pions as projectiles on heavy targets 
in the limiting fragmentation regime (Fiiab > 10 GeV). Of course, reverse 
kinematics is also possible in the latter case, as is the use of light ions with 
a corresponding energy deposition. While the use of light ions already leads 
to cavitation [2], the use of near symmetric projectile-target combinations 
leads to a significantly higher compression that is released into a high value 
of radial flow, which can cause the appearance of bubble-nuclei jSj. While 
the transient creation of bubble nuclei is of course interesting in its own right, 
it can also lead to a significant change in the observed fragment spectra [1]. 

We have introduced the percolation model of nuclear fragmentation al¬ 
ready two decades ago [SliaiTlIHl. It represents the nucleons by a collection 
of lattice sites on a three-dimensional lattice (usually simple cubic, but not 
necessarily) and their interactions by nearest-neighbor bonds. One can show 
[H] that this prescription approximately reproduces the Bethe-Weizsacker 
binding energy systematics. Energy is deposited via the random breaking 
of bonds, and sites still connected to each other are identified as fragments. 
The number of broken bonds is strictly proportional to the total energy 
deposited; this calculation can be done in a microcanonical way (fixed num¬ 
ber of bonds), a a canonical way, relating the bond breaking probability to 
the temperature via [HI Pb = 1 ~ 2r (3/2,0, i?/T) /t/T. We have used the 
Glauber approximation jTj or intranuclear cascades m or experimental data 
HU to estimate the energy deposition and residue size after pre-equilibrium 
emission as a function of impact parameter. 

The calculations on the right hand side of fig ^ (histograms) were per¬ 
formed with this percolation model and a subsequent multi-body Goulomb 
expansion. The observed fragment distributions are sensitive to the freeze- 
out density and thus enable us to determine this important parameter. The 
calculated fragment spectra are of Boltzmann-shape, but the slope of the 
spectra is not equal to the intrinsic temperature of the hot nuclear system 
at breakup HU- 

As a further sampling of the successes of the percolation model we show 
in fig. 121 on the left a comparison of the charge yields in the reaction of 
10.8 GeV pions with gold nuclei m and on the right an event-by-event 
reconstruction of the second moment of the fragment charge distribution as 
a function of multiplicity 113 HU for the inverse kinematics reactions of 1 


A GeV gold nuclei on a carbon target. 

The percolation model has a well-defined universality class in the infinite 
size limit. The cluster size distribution in the vicinity of the critical point is 
governed by the critical exponents a and r and follows the scaling law 

n{A,p) = aA~'^f{A^ ■ {p-Pc)) (for p « pc) (1) 

where A is the cluster size, p is the control parameter and pc its critical 
value. / is the scaling function with the property /(O) = 1. Thus, if one 
plots the log of the scaled yield, ns{p)s'^, versus the scaled control parameter, 
{p — Pc)s^, the cluster distributions produced in percolation calculations all 
fall on a straight line that crosses the point (0,1). 



Figure 2: Left: Charge yields from the reaction of high energy pions with gold, data 
from El; right: event-by-event distribution of the second moment of the fragment 
charge distribution as a function of the multiplicity including (upper curves) and 
excluding (lower curves) the largest fragment in each event, data from m 

One can also conduct an analysis of experimental fragment size distri¬ 
butions, using the value of the critical temperature and those of the critical 
exponents a and r as free parameters, then performing a optimization 
procedure to find the set of the values for which the scaling collapse of the 
data is optimized. This was done in HH for the 10.8 GeV pion-induced 
gold fragmentation data. The values of the critical parameters obtained are 
a = 0.5 ±0.1, T = 2.35 ±0.05, and = 8.3 ±0.2 MeV. However, we need to 
point out that this analysis is not entirely model independent, because one 

























needs to correct for the effects of sequential decays and feedings. Similar 
values for the critical exponents were reported in m- 

Recently, Zipf’s Law has been applied to the size of the largest cluster in 
nuclear fragmentation. Zipf’s Law in its original formulation applies to the 
frequency of the most frequent words in the English language. Zipf found 
m in the 1940s that if one ranks the words by their frequency, one finds that 
the frequency F has a power-law dependence on the rank r: F{r) = cr~^, 
with c = constant, A ~ 1). The same distribution was found much earlier 
by Pareto in ranking companies by their income EH- Similarly, ranking US 
cities by their population or citation analysis of scientific papers or many 
other distributions yields the same law. Recently, Watanabe m observed 
in numerical simulations that percolation clusters at the critical threshold 
also obey Zipf’s Law, and Ma nnum proposed to make use of this finding 
to detect a crossing of the critical threshold. 

One can derive m exact expressions for the average size of the largest 
clusters in a fragment distribution at the critical point, which is generated by 
a scaling theory that obeys eq. (©. Under the assumption that the system 
size V is large compared to the size of the emitted cluster, the probabilities 
Pk(i) to have k clusters of size A in a given event are Poissonian. And the 
probability that a cluster of size A is the largest in a given fragmentation 
event 

Pi (A) = [l-po{A)] ■po{>A) = [1 - • eI<(-4+A)-aC(r,i+y)]^ ( 2 ) 

with the normalization constant a = The probability that 

a cluster is the second largest is 

P2(A) = p> 2(A) •po(> A)-hp>i(A) •pi(> A) 

= [^ - Po{A) - pi[A)\ ■ po{> A) + [I - po(A)\ ■ pi{> A) (3) 

The probabilities for higher ranks follow from recursion relations. From this 
we obtain the exact expression for the average size of the largest clusters by 
summation over all A, (A^) = ^ ‘ TV (A). 

In fig. 01 we show the results (plot symbols) from these (exact!) simu¬ 
lations for different values of the critical exponent r and for a system size 
of 10,000 constituents. We take the ratio of the average size of the largest 
cluster to that of the cluster of rank r. If Zipf’s law were to hold, one would 
expect a strictly linear function with slope one, if plotting this ratio versus 
the rank (straight line without plot symbols). One can see that different 
values of r lead to different curves, and that none strictly adheres to Zipf’s 
law. 
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Figure 3: Average size of largest cluster over average size of cluster of rank r as a 
function of r for different values of the critical exponent r. The expectation from 
Zipf’s law is represented by the straight line without plot symbols. 


If one uses an expansion for our analytical results, one finds that the 
average cluster size as a function of rank r follows a more general Zipf- 
Mandelbrot distribution |22[ I23| . (A^) = c(r + k)~^, where the offset k is 
an additional constant that one has to introduce, and A is asymptotically 
approximated as a function of the critical exponent r, A(r) = l/(r — 1). 
These distributions are represented by the solid lines in fig. |31 

Thus the result of our calculations is that the cluster size distributions do 
not follow Zipf’s Law, but instead follow the more general Zipf-Mandelbrot 
distributions m- 
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